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f**") Aims. We constrain the past star formation histories of a sample of 25 distant (z ~ 0.7) luminous infrared galaxies (LIRGs) detected with the 
mid infrared cameras ISOCAM and MIPS onboard the ISO and Spitzer satellites. 

Methods. We use high resolution VLT-FORS2 spectroscopy in addition to a comprehensive library of 200,000 model optical spectra to de- 
rive Bayesian likelihood estimates of the star formation histories of these galaxies based on the analysis of Balmer absorption lines and the 
O 1 4000 A break. 

^—j , Results. The locus of distant LIRGs in the diagram defined by HSa and D4000 is roughly comparable to that of local LIRGs observed with 
(/} ■ IRAS, suggesting that no trend for an evolution is detected between the local and distant LIRGs. We obtain similar results when using either 

' the H8 or the H<5a Balmer absorption-line indices in combination with D4000. 
^ We compute a birthrate parameter (b=SFR/<SFR>) of 4±1, confirming that the distant LIRGs are currently experiencing a major phase of star 
V T ~j ■ formation. The most likely duration of the bursts is 0. lO+jJ^ Gyr, during which the LIRGs produce ~5-10% of their current stellar mass. No 
evidence is found for successive starbursts on the scale of a few times 10 7 yr, such as those predicted by some numerical simulations of major 
£2 . mergers. However, the large number density of those galaxies suggest that they could experience between two and four LIRG phases until the 
' present epoch. This scenario is not consistent with the formation of the z~0.7 LIRGs through the continuous star formation characterizing 
isolated spiral galaxies as has been argued independently based on their morphology. Instead, minor mergers, tidal interactions or gas accretion 
remain plausible triggering mechanisms for more than half of the distant LIRGs which do not harbor the morphology of major mergers. 
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1. Introduction 

Luminous infrared (IR) galaxies, i.e. galaxies radiating more 
than ~ 90 % of their light above 5 fim, have been suggested to 
provide important constraints on galaxy formation and evolu- 
tion. They are considered to be the main cause of the cosmic 
infrared background (ORB) and major contributors to the evo- 
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lution with redshift of the cosmic star formation rate (CSFR) of 
galaxies (Elbaz et al. 2002, Chary & Elbaz 2001, see also the 
review by Lagache et al. 2005). The role of luminous IR galax- 
ies (LIRGs for 10 n L Q < L m = L[8-1000jtm] < 10 12 L o )and 
ultra-luminous IR galaxies (ULIRGs for Lir > 10 12 L0) in the 
local universe could be neglected since they only produce about 
6 % of the integrated IR emission of local galaxies, hence 2% of 
their bolometric luminosity (Soifer & Neugebauer 1991). But 
they rapidly evolve with redshift and dominate the CSFR above 
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z ~ 0.6 (Le Floc'h et al. 2005). IRAS extragalactic surveys al- 
ready suggested that their number density evolved rapidly with 
redshift (~ (1 + z ) 7 - 6±3 - 2 up to z ~ 0.2; Kim & Sanders, 1998). 
This rapid evolution was later on confirmed up to z ~ 1 with 
ISOCAM and ISOPHOT onboard ISO at 15 and 170 /mi (see 
Elbaz et al. 2005 and references therein), above z ~ 2 using 
SCUBA at 850/on (Smail et al. 2001, Blain et al. 1999) and in 
the intermediate redshift range with MIPS onboard Spitzer at 
24 /mi (Chary et al. 2004, Papovitch et al. 2004, Le Floc'h et 
al. 2005). However, little is known about the characteristics of 
the starbursts themselves, e.g. the amount of stars born during 
the burst, the burst duration, and even less about the physical 
processes responsible for the intense activity of these galaxies. 
One of the major reasons for this lack of information comes 
from their nature itself, i.e. the strong dust obscuration of their 
optical light. 

In a previous paper (Liang et al. 2004, hereafter Paper I), 
we presented an analysis of the emission line properties of 
the galaxies. The star formation rates (SFR) derived from the 
Balmer emission lines (Ha and/or H/3, plus Hy to derive dust 
attenuations) were corrected for dust attenuation and found to 
be consistent with the ones derived from the mid IR (MIR) 
using the technique described in Chary & Elbaz (2001). This 
study showed that LIRGs in general are not completely ob- 
scured by dust and that the use of high-resolution optical spec- 
troscopy (AA/A = 2000, in the rest-frame of the objects) could 
be used to derive intrinsic luminosities, hence SFR, in rough 
agreement with the IR-derived SFR, by minimizing the con- 
tamination by sky emission lines and allowing to better cor- 
rect for underlying photospheric absorption lines. However, the 
consistent derivation of the signal-to-noise ratio on the intrin- 
sic luminosities lead to large uncertainties on the measured vi- 
sual attenuation. Moreover, there is evidence for some com- 
pletely obscured star formation as found in the most luminous 
objects studied in Flores et al. (2004) or in Cardiel et al. (2003). 
The limited statistics of those studies clearly calls for an ex- 
tention of the sample of distant LIRGs, with good S/N on the 
optical continuum and high spectral resolution, to robustly de- 
termine which fraction of the star formation taking place in 
LIRGs and ULIRGs is completely obscured by dust. However 
Hopkins et al. (2003) showed that SFR([OII]), SFR(1.4 GHz), 
and SFR(FIR) are in very good agreement for a larger sample 
of local infrared galaxies detected with IRAS and spectroscop- 
ically observed with the Sloan Digital Sky Survey (SDSS). 

In the present paper, we wish to address the problem from 
another angle: stars less massive, hence with longer lifetimes, 
than those responsible for the emission lines standardly used to 
derive the optical SFR of galaxies can escape their parent giant 
molecular cloud (GMC) and their spectral signature might be 
used to derive key parameters concerning the starburst. The Ha 
light used to derive a SFR is dominantly produced by the ioniz- 
ing photons arising from stars more massive than ~10 M Q , with 
lifetimes shorter than 3 Myr. Those stars never escape their par- 
ent GMC (average lifetime of 10 Myr) and the dense regions of 
very strong extinction, contrary to the A and F stars which are 
the main contributors to the Balmer absorption lines and the 
4000Abreak. 



We use these signatures of the optical continuum to com- 
pare distant LIRGs to nearby IRAS galaxies or synthetic spec- 
tra generated with the Bruzual & Chariot (2003) model. We 
find a signature of the starbursts in the D4000-H5 diagram and 
use it to derive the burst properties. 

Paper I indirectly confirmed the strong role played by 
LIRGs in the CSFR history with the derivation of gas metal- 
licities in distant LIRGs twice smaller than the one measured 
in present-day galaxies of equal absolute B band magnitude. 
This result suggested that these galaxies produced about half 
of their metals between z ~ 1 and today. It was also suggested 
that such a large metal production, as well as the large contri- 
bution of LIRGs to the CSFR and CIRB, could not result from 
a single burst phase in the galaxies harboring LIRG phases and 
hence that those galaxies must have experienced a series of 
LIRG phases in their lifetime. This possibility is tested in the 
present paper. 

Finally we note that, based on optical spectra, no evidence 
was found in Paper I for a dominant contribution from ac- 
tive galactic nuclei (AGNs) in the sample of 76 distant LIRGs 
for which a spectroscopic redshift was determined. This con- 
firmed the previous result from Fadda et al. (2001) that AGNs 
were contributing to less than 20 % of the MIR light of distant 
LIRGs, as shown by their soft to hard X-ray radiation measured 
with the Newton and Chandra X-ray observatories. 

In this paper, we study the stellar spectra of a sample of 25 
LIRGs with a median redshift of z = 0.7 using high-resolution 
spectroscopy with the FORS2 instrument at the VLT (AA/A = 
1200, equivalent to 2000 in the rest-frame of the galaxies). This 
sample is smaller than in Paper I because higher signal-to-noise 
ratios are required to study the continuum emission in compar- 
ison to the emission lines. All galaxies are detected at 15//m 
with ISOCAM and the 1 1 galaxies located in one of the three 
fields are also detected at 24 //m with the MIPS camera onboard 
Spitzer (Papovich et al. 2004, Elbaz et al. 2005). We show that 
both indicators provide consistent estimates of the total IR lu- 
minosity of the galaxies, hence also the SFR (see also Elbaz et 
al. 2005, Marcillac et al. 2005 for a more detailed study). 

Sect. 2 presents the sample selection and the wavelet de- 
composition technique that we used to analyze the spectra. 
Sect. 3 describes the method used in this paper to study the SFR 
history of the galaxies, namely the Balmer absorption line H6a 
(4101 A) versus 4000 A break position of the galaxies, as previ- 
ously done by Kauffmann et al. (2003) for the Sloan Digital Sky 
Survey (SDSS). We extended the method used in Kauffmann et 
al. (2003) to the high order Balmer line H8 (3889 A) and H9 
(3835 A)which are easier to detect in distant galaxies due to the 
k-correction. The comparison of local and distant LIRGs in this 
parameter space is discussed in Sect.4 while the model used to 
generate Monte Carlo realizations of 200,000 spectra with dif- 
ferent star formation histories is presented in Sect.5. The results 
are presented in Sect. 6 and discussed in the last Section. 

Throughout this paper, we will assume H„= 75 km s _1 
Mpc -1 , Q ma tter= 0.3 and D. A = 0.7. 
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Table 1. Description of the distant LIRGs sample. The technique used to compute the equivalent widths of H5, H8, H9, H10 and 
the 4000 A break (D4000) is described in Sect.El L 1R and SFR IR were derived using the Chary & Elbaz (2001) technique. 



2. Sample selection and data reduction 

2.1. Sample selection 

The sample of distant LIRGs was selected from three deep 
ISOCAM 15 /im surveys and is described in Paper I. 105 galax- 
ies were selected on the basis of their 15^m flux density in 
three different regions of the southern hemisphere, hence avoid- 
ing strong contamination by cosmic variance: the ISOCAM 
Ultra-Deep Surveys in the FIRBACK (UDSF) and ROSAT 
(UDSR) fields and the CFRS 3 h field. All three fields were se- 
lected for their low cirrus contamination and high galactic lati- 
tude. The UDSF (9'x9') is located at the center of the "far in- 
frared background" survey at 175 /mi with ISOPHOT onboard 
ISO ('FIRBACK'; Puget et al. 1999, Lagache & Dole 2001). 
The UDSR is centered at the position of a deep ROSAT survey 
(Zamorani et al. 1999). The UDSR and UDSF are both close to 
the position of the so-called 'Marano Field' originally selected 
for an optical survey of quasars (Marano et al. 1988), but are 
separated by 2 1 'with respect to each other. The third field is one 
of the Canada France Redshift Survey fields (CFRS 3 h ) com- 
bining deep infrared, optical and radio data as well as spectra 
from the MOS multiobject spectrograph on the 3.6m CFHT 

ISOCAM sources were selected to span the whole flux den- 
sity range of the three surveys whose 80 % completeness lim- 
its are 150yuJy for the UDSF and UDSR and 250 /Jy for the 
CFRS 3\ while the detection limits are 50/Jy (UDSF, UDSR) 



and 170/vJy (faintest sources in the regime below complete- 
ness). No optical selection was applied so that the magnitudes 
of the sources range from Rab- 19.4 to 23.7. At about the 
same epoch (November 2003) than the VLT-FORS2 observa- 
tions, a 24//m survey was being performed with the MIPS 
camera onboard Spitzer during the MIPS commissioning phase 
(IOC/SV). This survey covers the whole UDSF field and all 
eleven ISOCAM-15/mi sources were detected at 24fim, i.e. 
45 % of the sample. The scan map AOT was used, with an 
half-array overlap to cover about 1300 Sq. Arcmin with high 
redundancy (20) and to get an integration time per sky pixel of 
about 230s (Papovich et al. 2004). The data were reduced us- 
ing the Spitzer Science Center Pipeline and the BCD products 
(Basic Calibrated Data, Pipeline version S 10.0.3) were down- 
loaded from the Spitzer archive 1 . PSF-fitting photometry was 
performed using DAOPHOT (Stetson 1987) with IRAF 2 . 

From a total of 105 galaxies observed with FORS2 at the 
VLT (Paper I), we detected 3 stars, 13 galaxies were too faint 
for a redshift determination and within the remaining list, 1 3 
other galaxies were not detected with ISOCAM. The result- 
ing list of 76 galaxies with a MIR detection is divided into 34 



1 http://archive.spitzer.caltech.edu, PID: 718 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement 
with the National Science Foundation 
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Fig. 2. Comparison of the VLT-FORS2 raw spectrum (grey 
line) of the distant LIRG UDSR23 (z=0.7094, 51.8 M yr 1 ) 
to its wavelet decomposition resulting from the combination of 
the four wavelet scales 4-7, i.e. from 1 1 .2 to 89.6 A(black line). 
While the absorption lines are clearly distinguishable, the emis- 
sion lines such as [OIII],[OII],HyS,Hy are partly diluted by this 
technique. 



"normal" galaxies, 36 LIRGs and 6 ULIRGs, where "normal" 
galaxies consist in all galaxies except IR luminous galaxies (i.e 
LIRGs and ULIRGs). The median SFRs associated to each to- 
tal IR luminosity bin are 4, 54 and 196 M Q yr~' respectively. 
Because good quality spectra are required to study the stel- 
lar absorption lines (S/N> 3 on the continuum per resolution 
element), the final sample that we study in the present paper 
consists of 25 LIRGs with IR luminosities ranging between 
10 11 and 1O 12 L . The median redshift of the final sample is 
z=0.7. The measured properties of the galaxies are summarized 
inTableQ] 

The total (8-1000 pm) IR luminosities, Lir, were derived 
using the library of template SEDs built by Chary & Elbaz 
(2001), as in Elbaz et al. (2005). We also computed Lir from 
the Dale & Helou (2001) library following the technique de- 
scribed in Marcillac et al. (2005) and found a median value for 
Lir 10 % lower than with the previous library and with an rms 
of 17 %, hence both techniques provide consistent luminosities. 
We then compared the Lir derived with 15 pm and/or 24 fim 
flux densities for the 1 1 galaxies detected with both ISOCAM 
and MIPS. Both libraries of template SEDs provide consistent 
determinations of Lir using both measurements (with an rms 
of 30 %, see also Elbaz et al. 2005). The median Lir derived 
from MIPS is 10 % lower than the one derived from ISOCAM 
using both libraries which suggests a possible variation of the 
MIR spectra of LIRGs as a function of redshift (see Marcillac 
et al. 2005). 

2.2. Observations and data reduction 

The observations were performed during three nights with 
FORS2 on the ESO-VLT with the combination of the grisms 
R600 and 1600 (3 hours per grism) to cover the wavelength 
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Table 2. Definition of the wavelet scales used in the wavelet 
decomposition of the VLT-FORS2 spectra. The spectral reso- 
lution is equal to the highest resolution of o.7x2 waveletscale A. 
See description in the text for the construction of the different 
wavelet transforms at each scale. 



range 5000 to 9200 A at a resolution of 5 A (R=1200). At the 
median redshift of the objects of z ~0.7, the resolution is equiv- 
alent to 3 A (R=2000). 

Because absorption lines are wider than emission lines, it 
is possible to increase the signal-to-noise ratio on the absorp- 
tion lines by working at a lower frequency than the raw spec- 
trum. The physical origin of the width of the absorption lines 
results from the complex combination of the internal dynamics 
of the stars and the global dynamics of the galaxy. In order to 
optimize the line extraction, we therefore decomposed the raw 
spectra into eight different wavelet scales (Table |2) using the 
undecimated (keeping an identical sampling in each wavelet 
scale) wavelet transform (a trous algorithm; Starck & Murtagh 
1994, Starck, Pantin & Murtagh 2002). We optimized the S/N 
ratio on the absorption line features by selecting the best com- 
bination of wavelet scales. In the wavelet space, the first scale 
(highest frequency), that we index as scale 0, corresponds to 
features of size 0.7 A, while features at scale j have a size of 
0.7x2 J A. Note that the lowest frequency scale is equivalent to 
the raw spectrum smoothed by a B-spline (equivalent to a gaus- 
sian) of width 89.5 A, i.e. the baseline, while at higher frequen- 
cies, each wavelet scale T is equal to the difference between 
the raw spectrum smoothed at the scales T and 'i-1'. Hence 
the sum of all wavelet scales plus the baseline (here at scale 7) 
is exactly equal to the initial raw spectrum. 

In order to determine the number of scales to take into ac- 
count in the decomposition, we started from the lowest resolu- 
tion in the wavelet space, equivalent to a spectral resolution of 
89.5 A. We then kept doubling the frequency level as long as 
the S/N ratio was increased. This happened at the wavelet scale 
'4' (equivalent to 11.2 A), hence we used the co-addition of the 
4,5,6 and 7 wavelet scales (from 11.2 to 89.6 A, increasing the 
frequency by a factor 2 for each scale) in order to reconstruct 
spectra devoid of high resolution noise. The steps of the de- 
composition are shown in the Fig[2for a LIRG, UDSR23 (z= 
0.7094), located at the median redshift of the sample. The spec- 
tra resulting from the combination of these four scales present 
the advantage of seeing the same spectral resolution as the one 
used by Bica & Alloin (1986). The final result is compared to 
the raw spectrum in the Fig. 13 
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Fig. 1. Wavelet decomposition of the VLT-FORS2 spectrum of a distant LIRG (UDSR23) located at a redshift of z=0.7094) and 
forming stars at a rate of SFRir= 51.8 M yr~' . The raw spectrum (thick grey line) is compared to its wavelet decomposition (thin 
dark line), starting from scale 7, corresponding to the spectrum smoothed with a B-spline (equivalent to a gaussian) at the lowest 
resolution of 89.6 A, then adding up the wavelet scales by increasing the spectral resolution from 44.8 A(scale 6) to 0.7 A(scale 
1), by dividing it by a factor 2 for each scale. The narrow vertical grey zones were masked in the analysis because of the presence 
of sky emission lines. The vertical dotted lines indicate the minimum and maximum wavelength for each absorption line. In this 
zoom, the absorption lines clearly appear as much shallower than the emission line in OIL 



Note that the wavelet decomposition can dilute emission 
lines and that we have checked that absorption lines were not 
affected by a similar effect. Being wider, they are naturally less 
affected by this technique. However, we quantified this effect 
using some stellar spectra extracted from STELIB (Le Borgne 
et al 2003) which resolution is about the same as the rest-frame 
one for the distant LIRGs. We added a white noise to the stellar 
spectra to reach a signal to noise ratio of S/N ~ 3 on the con- 
tinuum, interpolated the spectra to reach a resolution of 0.7 A 
and applied to them the same wavelet decomposition as for the 
distant LIRGs. The equivalent widths of the Balmer absorption 
lines determined before and after the wavelet analysis differ by 
less than 4%. We included this weak difference in the equiva- 
lent width uncertainties. 



2.3. Measurement of the Balmer absorption line 
indices and 4000Abreak 

In the following, we will compare the equivalent widths mea- 
sured for H6 (4101 A) and for the high-order Balmer absorp- 
tion line H8 (3889 A). We also considered using the H9 (3835 
A) line that we therefore also define in this section, but it ap- 
peared to be a less reliable tracer of the past star formation 
activity of galaxies, as we discuss in in the next section. 

We used the definition of the H6 pseudo-equivalent width 
indice as defined in the Lick system (Worthey & Ottaviani 
1997)(see table [3Jl. Since no Lick indices have yet been de- 
fined for the high order Balmer absorption lines, we adapted 
the windows defined by Bica & Alloin (1986) to the Lick index 
method for H8 and H9 : for each line, the index continuum, the 
blue and red bandpasses for each pseudo-continuum are sum- 
marized in table[3] The two latter lines present the advantage of 
being located at a lower wavelength, to be accessible to higher 
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Name 


Blue and red bandpasses (A) 


Index continua(A) 




HS A 


4041.60-4079.75 4128.50-4161.00 


4083.50-4122.25 


10 


H9 


3810-3820 3855-3865 


3825-3845 


8 


H8 


3855-3865 3905-3915 


3870-3900 





Table 3. Definition of the pseudo-equivalent width indices for 
the H6 (H6 A in Worthey & Ottavianni 1997) and for the two 
high order Balmer absorption lines. For the two last lines, we 
followed the same principle as for Lick indices while using 
windows defined Bica & Alloin (1986). 
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redshifts in the observed optical range and to be less affected by 
the overlying nebular emission lines at the same wavelengths. 

The H5 absorption line is surrounded by iron absorption 
lines which affect both its red and blue pseudo continua and 
are responsible for the negative EW measured for this line after 
a few Gyears, as discussed in the next section. In order to avoid 
such pollution from neighbouring lines, the red continuum of 
H8 was choosen to minimize the sensitivity to the Call H (3933 
A)and K (3966 A) lines while the 3855-3865 A region is not 
known to be affected by metallic lines. As a result this index 
is weakly affected by metallicity. H9 is the least polluted line 
because no metallic line is located near it, but it is somewhat 
fainter than H8 and it is surrounded by the strong absorption 
of H8 and H10 which makes the two pseudo continua more 
difficult to define. We have primarily used H6a and 7/8, which 
are better determined for most of our spectra but also checked 
that we obtained consistent results based on the H9 and H6 
lines when they were available. Note that we will discuss sep- 
arately the results obtained with the Balmer absorption lines 
U6 A (4102 A) and H8 (3889 A). The advantage of this choice 
is that it provides two independent estimates of the parameters 
that we are deriving and can be used as a test of the robustness 
of the Bruzual & Chariot (2003) code that we are using. 

Note also that there is an ongoing debate about the possi- 
ble misinterpretation of the equivalent width of the H6a line 
because of metallicity ratios which could affect its neighboring 
regions, hence its associated pseudo-continua (see Thomas et 
al., 2004, Korn et al., 2005). On the one hand, HSa presents the 
advantage to allow the comparison with studies of local galax- 
ies such as SDSS galaxies (Kauffmann et al. 2003) while H8 is 
usually not available for local galaxies because it lies in a bluer 
region of the spectrum. On the other hand, the H8 line is not 
known to be affected by neighboring metallic lines and as it is 
located in a bluer region of the spectrum, it is easier to measure 
for distant galaxies (less polluted by sky emission lines). 

Before measuring these absorption features, we have cor- 
rected them from the overlying nebular emission line whenever 
possible as indicated in the paper. The nebular emission lines 
are not detected directly from the spectra, because they are too 
faint, but we computed their emission based on the observed 
Hp and H y emission lines (and H a for the low-z galaxies) as- 
suming a line ratio corresponding to a case B recombination 
for electron densities < 10 4 cm -3 and temperatures T e ~ 10 4 
(Osterbrock, 1989). The Balmer emission line ratio was also 
used to compute the attenuation of these lines before subtract- 
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Fig. 3. Time evolution of the H8, H9 and H6a indices and 
4000Abreak (D4000) for a single stellar population synthe- 
sized with the GALAXEV code (Bruzual & Chariot 2003). 
Four different metallicities were used in each plot demonstrat- 
ing the marginal effect of metallicity on the three indices and 
on D4000 before 1 Gyr: Z= 0.004 (20 % solar, dashed line), Z= 
0.008 (40 % solar, dash-dotted line), Z= 0.02 (solar, solid line), 
Z= 0.05 (2.5 times solar, dotted line). From top to bottom and 
left to right: a) H8 evolution, b) H9 evolution, c) HSa evolution, 
d) D4000 evolution. 

ing them from the absorption lines measurements. The compu- 
tation and values for these attenuations can be found in Paper 
I. 

Following the same strategy as Kauffmann et al. (2003), we 
used the method explained in Bruzual (1983) and the blue and 
red bandpass definitions introduced by Balogh et al. (1999). 
The latter are narrower than the ones originally defined by 
Bruzual (1983) and present the advantage of being less sen- 
sitive to reddening effects. 

H8 and D4000 values for each distant LIRG are summa- 
rized in Tabled The scientific interpretation is discussed in the 
following section. 

3. 4000 A break and Balmer absorption lines as 
tracers of the recent star formation history 

In this section, we describe how Balmer absorption lines and 
the 4000 A break trace the recent star formation history of 
galaxies. For this purpose, we have synthesized a single stellar 
population (SSP) using the latest version of the 'GALAXEV 
code from Bruzual & Chariot (2003). For the moment, we 
do not include dust attenuation to simplify the discussion. 
Note however that the wavelength range over which the equiv- 
alent widths and the 4000 A break are measured is small 
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which implies a marginal correction due to dust attenuation. 
Nonetheless, we will include dust attenuation in the Monte 
Carlo realizations that we will compare to the LIRGs and dis- 
cuss its effects on our findings. 

This version of the GALAXEV code includes the spectral 
library STELIB (Le Borgne et al. 2003) whose spectral reso- 
lution is 3 A from 3200 to 9500 A, which is comparable to the 
present spectra in the rest-frame of the galaxies. Fig.|3]presents 
the evolution with time of the high order Balmer absorption 
lines H8 and H9 as well as the H5a absorption line and the 
4000 A break. The four lines in each plot correspond to four 
different metallicities with the following metal mass fractions 
(total mass in elements heavier than hydrogen and helium over 
the mass in hydrogen): Z= 0.004 (20 % solar, dashed line), Z= 
0.008 (40 % solar, dash-dotted line), Z= 0.02 (solar, solid line), 
Z= 0.05 (2.5 times solar, dotted line). 

The D4000 is the sudden onset of stellar photospheric opac- 
ity shortward of 4000A. It reflects the mean temperature of the 
stars responsible for the continuum: the metals located in the 
atmosphere of O and B stars are more ionized and produce a 
weaker opacity, hence a smaller 4000 A break, than those in 
cooler stars (Bruzual 1983, Poggianti & Barbara 1997, Gorgas 
et al 1999, Kauffmann et al 2003). Hence D4000 keeps in- 
creasing as a function of the aging of the stellar population 
(see Fig. |3Jl). D4000 is sensitive to metallicity as it is shown 
in tab. [3jl where it varies by more than 20 % after a few billion 
years, or when it is larger than 1.6. As we will see in the next 
section, the distant LIRGs have D4000= 1.2+0.07, and there- 
fore for them metallicity effects are negligible. Even if the slope 
of D4000 versus age is less flat than for young population, i.e. 
7 Gyears, it is flat enough to provide uncertain stellar ages if 
used alone. 

In order to trace back the recent star formation history of 
galaxies, it is therefore necessary to use another tracer of stellar 
age such as the Balmer absorption lines, which exhibit a steep 
slope as a function of stellar ages in this range of ages (see 
Figs.|5t,b,c). Balmer absorption lines are dominantly produced 
by the atmosphere of A to F stars. However O and B stars, 
which do not exhibit strong absorption lines, indirectly affect 
them by increasing the continuum level and therefore diluting 
them, which explains the flat values for the equivalent widths 
of H8, H9 and H6a in the first few millions years (lifetime of 
O and B stars). The rapid increase that follows is produced by 
the dominant role of A and F stars, which then disappears after 
~ 0.5 Gyr producing the rapid decline of the equivalent widths 
in Fig. |5Ji, b and c. Here again, it is worthwhile noticing the 
marginal role played by the metallicity on the evolution of the 
Balmer lines EW with time. The EW varies by less than 20 % 
a few billion years after the burst as a function of metallicity. 

4. Comparison of local and distant LIRGs 

A sample of 401 local (z <0.25) LIRGs detected with IRAS, 
with optical spectra from the SDSS and emission lines typical 
of star forming galaxies, as opposed to AGNs, was identifed by 
Pasquali et al. (2005). The locus of the D4000 and EW (U6 A ) 
measured for the distant LIRGs (green filled circles) is com- 
pared to that of local LIRGs (blue triangles) and of the field 
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Table 4. Median and 68 % dispersion (around the median) of 
the D4000 and H6a indices for local and distant galaxies. 



SDSS galaxies (black dots) in Fig. |4] First, note the concen- 
tration of both LIRG populations at low D4000 and high H6a, 
which suggests that local and distant LIRGs share a similar 
recent star formation history. We could not produce similar fig- 
ures for high order Balmer absorption lines for this compari- 
son since they are too blue to be accessible in the SDSS spec- 
tra. Their position in the D4000-H(5a diagram indicates that the 
young stellar population that is producing the large IR luminos- 
ity is not completely obscured by dust in the optical because 
the median D4000 value of 1 .2 for these galaxies correspond to 
stellar ages lower than 1 Gyr, which is much below the ages of 
these galaxies (see Table|4}. This already suggests that even in 
these dusty galaxies the optical spectral signatures can be used 
as a tracer of the recent star formation history. The quasi ab- 
sence of distant LIRGs above D4000= 1.25, where half of the 
local LIRGs lie, suggests that the distant LIRGs are younger 
than the local ones. 

The relative proportion of young and old stellar popula- 
tions can be studied in these galaxies by comparing the stellar 
masses, derived mostly from the old stellar population dom- 
inating in the near infrared range, with their D4000 values. 
Massive galaxies are generally older as shown by Fig.|5j where 
the stellar mass of SDSS galaxies (black dots) increases with 
D4000. The position of the local LIRGs (orange filled dots) in 
Fig. |5] suggests that they are massive galaxies which were lo- 
cated on the right side of the plot before the burst and which 
were shifted to the left during the burst, which decreased their 
D4000 value. This figure reinforces our interpretation that the 
low value of D4000 for the LIRGs is due to the addition of 
a young population on top of an older population. The com- 
parison of field SDSS galaxies with IRAS LIRGs in four bins 
of D4000 from 1 to 1.4 (bin size 0.1, 68 % error bars) shows 
that local LIRGs exhibit systematically lower D4000 values for 
their stellar masses. 

The incompleteness of the photometric data for our sample 
of galaxies prevented the determination of their stellar masses. 
However, another sample of mid IR selected LIRGs at z — 
0.7 located in the Hubble Deep Field South was studied by 
Franceschini et al. (2003) who computed their stellar masses 
using a Salpeter IMF and a combination of single stellar popu- 
lations with ten different ages, to fit their UV-optical-NIR spec- 
tra of IR luminous galaxies. From their Table 6, a total of 14 
LIRGs possess a spectroscopic redshift between z=0.4 and 1 .2, 
and 7 more have a photometric redshift in this range. After con- 
verting Franceschini's values to H = 75 km s _I Mpc -1 , we find 
a median stellar mass of 7xl0 ll) M for the 14 galaxies with 
a spectroscopic redshift. Including the less robust photomet- 
ric redhistfs only changes this value to 6x10 10 M Q . Combined 
with a median D4000~ 1 .2, this places the distant LIRGs in a 




Fig. 4. H6a versus D4000 diagram for the SDSS galaxies (black points) and distant LIRGs of the present sample (green filled 
circles with error bars). The sub-population of SDSS galaxies detected with IRAS are marked with triangles and separated into 
galaxies with optical spectra with (empty red triangles) and without (filled blue triangles) an AGN signature. 



very similar locus as the local LIRGs (see Fig. [5}- Their low 
D4000 values are therefore also representative of the young 
stellar population of the burst superimposed on top of an older 
stellar population dominating the stellar mass of the galaxies. 



5. Description of the model 

We used the Bruzual & Chariot (2003) model to synthesize a 
series of 200,000 Monte Carlo realizations including various 
star formation histories. Although the technique is similar to 
the one used in Kauffmann et al. (2003), it presents two differ- 
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D4000 

Fig. 5. Stellar mass as a function of D4000 for the SDSS galaxies (black points) and distant LIRGs of the present sample (large 
filled blue circle). The median stellar mass for the distant LIRGs was derived from Franceschini et al. (2003) as discussed in 
the text. The sub-population of SDSS galaxies detected with IRAS and without an AGN signature in their optical spectra (star 
forming galaxies) are marked with orange filled circles. The median and 1 -<x error bars for the field SDSS galaxies are represented 
with thick light-blue crosses, while the thin dark-blue crosses are for the sub-sample of star forming SDSS-IRAS galaxies. 
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Table 5. Description of the range of values used as priors for 
the simulations. 



ences related to the populations of galaxies that we are studying 
here, i.e. dusty starbursts. The priors of the model were set to 
include a range of dust attenuations and to include a larger frac- 



tion of starbursting galaxies. We will discuss the effects of both 
modifications in the following. 

Each star formation history was modeled with eight param- 
eters: 

- the age of the galaxy since its formation, tf orm . We im- 
posed a maximum age of 7 Gyr, which is the age of the universe 
at the median redshift of the galaxies, z ~ 0.7. 

- the parameter y for the underlying "continuous" SFR (as 
opposed to the starburst itself): SFR cont (t) oc e _7t 

- the SFR of the burst itself, SFR burst , is not a parameter 
but results from the combination of two parameters: the burst 
duration, tb, and the stellar mass fraction produced during the 
burst, /b- Note that fa computed in the model is the ratio of 
the stellar mass formed during the burst, M bulst , over the total 
stellar mass formed through continuous star formation during 
tform, M™ nt . Here M^ ont is not corrected for the mass returned to 
the interstellar medium by evolved stars and is therefore larger 
than the actual stellar mass of the galaxy, M+. In the follow- 
ing, we will call, f^ ff , the "effective burst stellar mass fraction", 
computed with M* instead of M™ nt . 
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D4000 D4000 D4000 

Fig. 6. Influence of the attenuation, velocity dispersion and metallicity for a typical star formation history synthesized with the 
GALAXEV code (Bruzual & Chariot, 2003). In all three plots, the loops correspond to a burst of star formation which started 
3 Gyr (big black circles) after the galaxy's birth, a) Influence of the attenuation r v : for a continuous star formation ( y=1.0, 
Vdisp=200km.s~', yi/=0.3, Z=Z Q ) and a single burst ( tb=0.1 Gyear and fb=10 % ); the solid line : t v =1.0. The dotted line : 
r v =2.0. The dashed line : t v =3.0. The long dashed line : r v =4.0. b) Influence of the velocity dispersion Vdisp : for a continuous 
star formation ( y=1.0, t v =3.0, /i=0.3, Z=Z q ) and a single burst (tb=0.1 Gyear and fb=10%). The solid line : Vdisp=100km.s~'. 
The dotted line : Vdisp=200km.s~'. The dashed line : Vdisp=300 Km.* -1 , confluence of the metallicity Z : for a continuous 
star formation ( y=1.0, Vdisp=200km.s~ 1 , r v =3.0 and fi=0.3 ) and a single burst (tb=0.1 Gyear and fb=10 %). The solid line : 
Z=O.5Z . The dotted line : Z=Z . The dashed line : Z=2.5Z0. 



- the dust attenuation parameters, and ry, defined as in 
Chariot & Fall (2000), i.e. a differential attenuation function 
of the age of the stars (see Eq. ^and Fig. [(J). Stars younger 
than 10 7 years, the typical lifetime of a giant molecular cloud 
(GMC), are assumed to be embedded in their dusty parent 
GMC with an optical depth ry. Older stars are supposed to 
have escaped their parent GMC and their ambient medium's 
optical depth is assumed to be ju times smaller. This attenua- 
tion law was chosen because it reproduces well the correlation 
observed for UV selected starbursts of the FIR over UV ratio 
with the UV slope, ft, as well as the Lh„ over Lhs ratio. The 
optical depth for stars younger than 0. 1 Gyears, Ty, was set be- 
tween 1 and 4 based on the observed values for LIRGs (Liang 
et al. 2004, Flores et al, 2004), while /i was allowed to vary 
from 0.1 to 1. 

^observed _ ^intrinsic ^ e T[ 

t a = t v x (5500M) - 7 u < 10 7 yr (1) 

t a = n x t v x (5500/A) - 7 U > 10 7 yr 

- the velocity dispersion, V^p. We allowed a large range of 
velocity dispersions to test the effect of a broadening of the H6a 
absorption lines that would mimick intrinsically large width of 
the absorption lines due to the combination of the internal dy- 
namics of the stars with the mass of the galaxy. The influence 
of the velocity dispersion on the position of a galaxy in the 
H6a-D4QQ0 diagram was found to be negligible (see Fig.^J)). 

- the metallicity, Z. The influence of the stellar metallic- 
ity is stronger than that of the velocity dispersion but remains 
smaller than the error bars on the measured H5 and D4000 (see 
Fig.0;). 

Table [5] summarizes the range over which those eight pa- 
rameters were chosen, the so-called priors, which represent the 
key adjustment when using a Bayesian statistics. The goal of 



these simulations is not to constrain all eight parameters and 
indeed some of them produce competitive effects or cannot be 
disentangled: 

- a galaxy with a given age and y parameter for its SFR will 
be located at the same position in the H8-D4000 diagram than 
a younger galaxy with a larger value for y. 

- the two parameters used in the attenuation law, ju and tv- 
Fig- El shows that galaxies are shifted towards the upper-right 
of the H<5^-D4000 (or equivalently H8-D4000) diagram with 
increasing dust attenuation. This is an important point to note 
when comparing dusty galaxies to field galaxies: even though 
LIRGs overlap with a fraction of SDSS field galaxies in the 
Fig- E] their recent star formation history is interpreted differ- 
ently if their dust attenuation is taken into account. 

Instead, we will use a bayesian approach such as the one de- 
tailed in Kauffmann et al. (2003) to derive the following char- 
acteristics of the starbursts: 

1. the Scalo parameter, i.e. the ratio of present to past- 
averaged star formation rate. This parameter defines the rel- 
ative intensity of the present starburst. 

2. the burst duration, tb • 

3. the last burst stellar mass fraction, /b (and /g ff ). 

The bayesian statistics approach consists in determining a 
probability distribution function (PDF) for any one of these 
parameters, which consists in an histogram of the number of 
Monte Carlo realizations weighted by the probability function 
exp(~x 2 /2), where x 1 is defined in Eq.f2] The convergence of 
the technique is reached if the PDF is peaked and its shape pro- 
vides the precision of the determination. 

_ / H5 obs - H6 sim \ 2 I D4000 obs - D4000 sira \ 2 
Xl ~ \ cr(H6 oh *) J + \ cr(D4000 obs ) j <_) 
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6. Results 

6.1. Signature of the presence of a starburst in the 
distant LIRGs in the H8-D4000 diagram 

We generated three series of simulations with the same pri- 
ors but varying the percentage of galaxies experiencing a star- 
burst at the age of the simulated spectrum. In a first simulation, 
SIM1, only 20% of the Monte Carlo realizations include on- 
going starburst. The positions in the H8-D4000 diagram of the 
200,000 Monte Carlo realizations of SIM1 are illustrated with 
small dots in Fig.0 The light blue points on the lower-left part 
of the figure illustrate the position of galaxies with an ongoing 
burst of star formation. The upper dark blue points correspond 
to post-starbursts, i.e. galaxies having experienced a recent star- 
burst which ended less than 2 Gyr ago, and dark points to galax- 
ies with continuous star formation or at least no burst during 
the past 2 Gyears. The lines illustrate the tracks followed by an 
individual galaxy during its lifetime: a galaxy with purely con- 
tinuous star formation will follow the bold red line, where the 
numbers mark its age in Gyr (y=2, /i=0.3, Tv=3.0,Vdi sp =200 
km.i -1 ). If this galaxy was to experiment a burst of star for- 
mation when it is 3 Gyrs old, then it would follow one of the 
loops in light orange towards the left of the diagram for a total 
duration of 2 Gyrs and then go back to the track of a galaxy 
with pure continuous star formation, having lost the memory 
of its past starburst. During the loop, D4000 is at first the most 
affected by the presence of young stars, moving towards the 
left of the diagram, then it is H8 which increases moving up- 
wards in the diagram, before the galaxy slowly goes back to 
the continuous regime. The first part of the loop is marked 
with small lines indicating that the galaxy is still bursting, 
while the remaining part of the loop corresponds to the post- 
starburst regime. We have plotted a loop starting at an age of 2 
Gyears (light orange) and another starting at 3 Gyears (dark or- 
ange). In each case we computed the track of Tg=5xl0 7 years 
(dashed lines) and 10 8 years (plain line) starburst with a fraction 
of/ B = 15%. 

75 % of the distant LIRGs are located within the solid rect- 
angle in Fig. hence at the location of galaxies experiencing 
a starburst in the simulation. The remaining 25 % of the distant 
LIRGs lie in the dashed rectangle which includes the region of 
post-starbursts but also galaxies with successive or longer star- 
bursts. Hence, the results of the simulation from SIM1 simula- 
tion show that the observed galaxies are mainly coherent with 
being starbursting galaxies. 

As a second step, we designed a new simulation, with 80 % 
of the galaxies experiencing an ongoing starburst (SIM2), to 
better sample the locus of the distant LIRGs and therefore bet- 
ter study the properties of the starbursts themselves. A third 
and last simulation (SIM3) was generated to quantify the prob- 
ability that the distant LIRGs experienced a previous starburst 
during the last 2 Gyr. Half of the galaxies in SIM3, i.e. 100,000 
galaxies, have experienced a previous starburst during the last 
2 Gyr, which ended before the onset of the ongoing starburst. 

The individual positions of the distant LIRGs are compared 
to the SIM2 realizations in Fig. |U and Fig|9]in the H8-D4000 
and H(5 J 4-D4000 diagrams respectively. UDSR09, which is ly- 
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Fig. 10. Example of PDF obtained for 1 galaxy : a) The Scalo 
parameter (subsection 6.2). b Burst duration (subsection 6.2). 
c) Burst fraction (subsection 6.2). d) Effective burst fraction 
(subsection 6.2). 



ing at the bottom of the simulation is a clear outlyer. Liang et 
al (2005) studied carrefully this object and its optical spectrum 
shows strong metal absorption lines such as Ca H K, G-band, 
Mg H, Na D lines but weak Balmer absorption lines. The X-ray 
emission of this object, in addition is mainly associated to an 
AGN. As a consequence, an AGN can partially contribute to the 
MIR luminosity. Note that the observed galaxies are not iden- 
tical in both figures neither in numbers nor in identity because 
both H8 and H6a cannot be both measured for all individual 
galaxies. The dots correspond to the same simulated galaxies 
in both figures and in both cases. Note that the observed LIRGs 
do lie below the continuous star formation regime in both fig- 
ures, although in slightly different locations. These differences 
will be discussed in the subsection I6.5l when we will present 
the resulting PDFs. Note also that the distant LIRGs are dis- 
tributed in two populations in Fig. [8] one located at the bottom 
of the diagram and a second at the upper-left. The second pop- 
ulation corresponds to simulated galaxies which experienced a 
succession of two starbursts during the last 2 Gyrs. Such histo- 
ries are expected in the framework of major mergers of spiral 
galaxies with several encounters between two galaxies. 

6.2. Determination of the Scalo parameter, 
SFR/<SFR> 

The PDF obtained for the Scalo parameter of a typical galaxy 
with SFR/<SFR>~ 2.7+™ (1-cr) is represented in Fig. ITUh, 
In about 60% of the distant LIRGs, the PDF converge to- 
wards a determination of the Scalo parameter. For 10 out of 
17 galaxies for H5a (see Table|6}, we compute a median value 
of SFR/<SFR>= 3.0^ (1-cr). For 12 out of 22 for H8 (see 
Table Q we find SFR/<SFR>= 4+| (1-cr). These two results 
are consistent with each other, as illustrated by the Figs. II lh 
and c. The fact that the H8-D4000 diagram provides the tight- 
est constraint on the Scalo parameter is due to the fact that the 
sky background is lower in the bluer wavelength range of H8 
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Fig. 7. Location of the 200,000 Monte Carlo realizations using the parameters of SIM1 (see table |5} in the H8 versus 
4000A(D4000) break diagram. Light blue points (lower-left): starbursting galaxies. Dark blue points (upper part): post- 
starburst galaxies (galaxies having experienced a recent starburst which ended less than 2 Gyr ago). Black points: galaxies 
whith a continuous star formation in the two past Gyears. Bold red line: track followed by an individual galaxy with continuous 
star formation (numbers in squares= age in Gyr) generated with GALAXEV (y=2., /i=0.3, Ty=3.0,Vdi sp =200 km.s -1 ). Orange 
loops: effect of starbursts of 5xl0 7 (dashed line) and 10 8 (plain line) years starting after 2 Gyr (light orange) and 3 Gyr (dark 
orange) of continuous star formation. The line is hatched during the starburst phase. The size of the orange dots is proportional 
to the time counted in 10 8 years units after the beginning of the burst. The dark rectangle with dashed contours: area where 
the distant LIRGS are located. The dark rectangle with a solid line: area where are lying 75 % of the sample (see Fig.[S]for 
the position of the individual galaxies with error bars). Triangles: galaxies corrected for the underlying nebular emission line in 
H8 as derived from the observed lower order Balmer lines corrected for attenuation from the Balmer lines ratio. Stars: galaxies 
with negligible emission Balmer line in emission, hence uncorrected for H8 emission. Squares: corrected for the underlying 
H8 emission line using only one Balmer emission line available with dust attenuation computed from the IR (from the ratio of 
SFR(IR) over SFR(H/3)). Circles: no correction for the underlying potential emission line because no information is available on 
the Balmer emission lines. 



than in W5a- We will therefore use that range of values for the 
Scalo parameter in the next set of computations. The fact that 
the LIRGs are found to produce stars at a rate about four times 
larger than their averaged past SFR confirms that they are expe- 
riencing a major phase of star formation in their lifetime. This 
result is consistent with the large Lir and SFR(IR). 

As noted above, the H8-D4000 diagram provides the tight- 
est constraint on the Scalo parameter with SFR/<SFR>= 4^ °. 
Combined with the ongoing SFR measured from the MIR 
emission for the distant LIRGs (quoted in Table ^ f° r indi- 



vidual galaxies and with a median value of: SFR = 52 ^ M 
yr~' , the Scalo parameter allow us to compute the averaged 
past SFR of the distant LIRGs : <SFR>~ 13+3 M yr _1 . 



If we assume that the progenitors of the distant LIRGs 
formed stars at a constant rate equal to the averaged past SFR, a 
median stellar mass of 7xl0 1() M (see Sect.|3 is assembled in 
about 5.5+1.5 Gyr (light- weighted age), implying that the first 
dominant stellar populations formed at a redshift z ~ 4.5±1.5. 
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Fig. 8. Same as in Fig.[7]for the SIM2 Monte Carlo simulation of 200 000 model galaxies. The distant LIRGs are represented 
with green points with error bars. 
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Table 6. Results concerning the burst duration, the Scalo ratio SFR/<SFR>, and the burst fraction (/b) when using the H6a data. 
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Fig. 9. Location of the 200,000 Monte Carlo realizations using the parameters of SIM2 (see table |3J in the H6a versus 
4000A(D4000) break diagram. Light blue points (lower-left): starbursting galaxies. Dark blue points (upper part): post- 
starburst galaxies (galaxies having experienced a recent starburst which ended less than 2 Gyr ago). Black points: galaxies 
whith a continuous star formation in the two past Gyears. Bold red line: track followed by an individual galaxy with continuous 
star formation (numbers in squares= age in Gyr) generated with GALAXEV (y=4., fi=0.3, Tv=3.0,Vdi sp =200 km.i -1 ). Orange 
loops: effect of starbursts of 5xl0 7 (dashed line) and 10 8 (plain line) years starting after 2 Gyr (light orange) and 3 Gyr (dark 
orange) of continuous star formation. The line is hatched during the starburst phase. The size of the orange dots is proportional 
to the time counted in 10 8 years units after the beginning of the burst. Triangles: galaxies corrected for the underlying nebular 
emission line in H6a as derived from the observed lower order Balmer lines corrected for attenuation from the Balmer lines ratio. 
Stars: galaxies with negligible emission Balmer line in emission, hence uncorrected for HSa emission. Squares: corrected for the 
underlying H6a emission line using only one Balmer emission line available with dust attenuation computed from the IR (from 
the ratio of SFR(IR) over SFR(HjS)). Circles: no correction for the underlying potential emission line because no information is 
available on the Balmer emission lines. 



6.3. Determination of the burst duration 



6.4. Determination of the burst stellar mass fraction 



In the Hi5a-D4000 diagram, the burst duration, tb, can be con- 
strained for 1 1 galaxies (Table|6) for which we obtain a median 
of T B = 1.3!^ x 10 8 years, e.g. Fig.[l0p. Using H8, the tech- 
nique converges for 14 galaxies to tb = 1.0*J"g x 10 8 years. 
Here again, as for the Scalo parameter in the previous section, 
we find a consistent result using both H8 and H6a- Both indi- 
cate that the bursts are short-lived over about 0. 1 Gyr. The large 
dispersion between the various galaxies may result from differ- 
ent initial gas mass fractions or triggering mechanisms for the 
burst. 



The simulations converge to a value of /b for only 6/22 galax- 
ies in the H8-D4000 diagram and 10/17 in the H(5 A -D4000 (see 
the example of UDSR 23 in Fig. Hclb), The difficulty to deter- 
mine the burst stellar mass fraction is not surprising. It results 
from the fact that as soon as the young stellar population starts 
dominating the spectrum of a galaxy, it becomes nearly impos- 
sible to measure with good precision the fraction of the light 
due to the underlying older population. Indeed, the PDFs con- 
verge only for galaxies for which the fraction of young stars 
remains low (see Tables l6T7Y 
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Table 7. Results concerning the burst duration, the Scalo ratio SFR/<SFR>, and the burst fraction (/b) when using the H8 data. 
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Fig. 11. Distribution of the Scalo parameter and the burst dura- 
tion, a) The Scalo parameter in the H8-D4000 diagram, b) The 
burst duration in the H8-D4000 diagram, c) The Scalo param- 
eter in the H#a-D4000 diagram, d) The burst duration in the 
H(5a-D4000 diagram. 



Fig. 12. The lines are the same as in Fig. and in Fig. [8] The 
black circle are proportional to £20 defined before, while the 
grey one to £50. 

molecular gas observed in local LIRGs and ULIRGs (Sanders 
& Mirabel 1996). 



However, we can still derive the effective burst stellar mass 
fraction, /| ff , which is equal to the mass of stars produced dur- 
ing the burst, i.e. tbX SFR, divided by the total stellar mass 
of the galaxy, i.e. f form x <SFR>. Since r B = 1.0^| x 10 s 
years and SFR=52 *|f M yr , we obtain an effective burst 
stellar mass fraction of = 10+9 %. Note that a burst pro- 
ducing 10 % of a 7xl0 10 M galaxy, will convert 7xl0 9 M Q of 
molecular gas into stars, which is consistent with the mass of 



6.5. Testing the "multiple burst" scenario 

A sub-sample of the distant LIRGs lie on the upper-left part of 
the H8 and H6 A versus D4000 diagrams (Figs. HOIS- This re- 
gion is populated by Monte Carlo realizations of galaxies for 
which the ongoing burst of star formation was superimposed 
on a previous burst which ended less than 2 Gyr ago. In or- 
der to quantify the probability of such an occurence, we gener- 
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Fig. 13. The lines are the same as in Fig. and in Fig. [8] The 
black circle are proportional to fio defined before, while the 
grey one to t$o. 
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Table 8. Values of f2o (black circle) and tso (grey circle) ob- 
tained for the 22 galaxies in the H8-D4000 diagram, to quantify 
the probability of two successive starbursts (see text). 



ated a third simulation, SIM3, in which half of the realizations 
experienced two successive bursts during the last 2 Gyr. We 
then computed for each galaxy, using the same bayesian statis- 
tics, a probability for various possible durations between the 
two bursts. Two representative numbers are assigned for each 
galaxy in Tables[8]and|9] tio and tso- These values correspond to 
lookback times associated with a 20 and 50 % chance of finding 
a previous burst which ended ?2o and tso Gyr before the onset 
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Table 9. Values of f2o (black circle) and tso (g re y circle) ob- 
tained for the 17 galaxies in the H<5a-D4000 diagram, to quan- 
tify the probability of two successive starbursts (see text). 



of the ongoing one. For example, in the case of UDSR20, we 
obtain from the H8-D4000 diagram (see TableE) f 2 o = 0.4 Gyr 
and fso = 1 Gyr, which implies that there is a 50 % chance that a 
previous starburst occured 1 Gyr before the onset of the present 
one and a 20 % chance that the delay was only a 0.4 Gyr. Hence 
UDSR20 is a good candidate for two successive bursts. We did 
not consider higher probabilities or longer timescales because 
of the limited constraints that we can set on those parameters 
and because after about 1.5 Gyr, the memory of the previous 
burst is lost with this technique. 

Local LIRGs and ULIRGs are known to be predominantly 
triggered by major mergers (Borne et al. 1999, Sanders & 
Mirabel 1996) and numerical simulations of such mergers pre- 
dict that the two galaxies cross each other several times, po- 
tentially inducing a series of bursts separated by a few tens of 
million years (Mihos & Hernquist 1996). However, in distant 
LIRGs, the probability that a previous starburst occured less 
than 0.5 Gyr ago is nearly always lower than 20%. This re- 
sult must be considered together with recent evaluations of the 
morphological properties of distants LIRGs which also suggest 
that most of them are not produced in major mergers (Bell et 
al. 2005, Zheng et al. 2004, Elbaz & Cesarsky 2004). Bell et 
al. (2005) suggested that distant LIRGs could either be non- 
triggered phases in isolated spirals with larger gas masses, pos- 
sibly experiencing some infall, or minor mergers, where the 
dwarf galaxy responsible is not detected. Our determination of 
a burst duration of 0. 1 Gyr and a Scalo parameter of 4, seems 
to rule out the possibility that distant LIRGs are isolated spi- 
rals forming stars at a constant rate over a long duration. The 
starbursts may instead be triggered by tidal effects and minor 
mergers in regions of the universe where the local density of 
galaxies is enhanced as suggested by Elbaz & Cesarsky (2003) 
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or by the infall of intergalactic gas (Combes, 2005). Further 
kinematical studies of distant galaxies (see Flores et al, 2006; 
Puech et al, 2006) will help to distinguish between the various 
scenarios (mergers, gas infall) discussed here." 

7. Discussion and conclusions 

We have analyzed the star formation history of a sample of 25 
distant LIRGs (z = 0.7) that we derived from their stellar spec- 
tra (Balmer absorption lines and 4000 A break). The high order 
Balmer absorption line H8 and the H6a line provide consistent 
results although H8 is more adapted to distant galaxies being 
located in a bluer region of the spectrum, less affected by sky 
emission lines and which is measured for more distant objects 
than H^. Variations at the 30% level between the burst pa- 
rameters obtained using one or the other indicator suggest that 
some effects, such as abundance ratios which might affect the 
pseudo-continua surrounding the HSa line, should be taken into 
account to improve the models (see Thomas, Maraston & Korn 
2004, Korn, Maraston & Thomas 2005). 

The comparison of distant LIRGs, selected from ISOCAM 
and MIPS, onboard ISO and Spitzer, to local LIRGs, selected 
from IRAS and the SDSS, shows that both populations present 
similar spectral features and therefore suggests that they are 
experiencing comparable starburst phases. The fact that half of 
the local LIRGs present D4000 values larger than the maxi- 
mum D4000 of distant LIRGs indicates that the dominant non 
bursting stellar population is younger for distant LIRGs, as ex- 
pected. 

The first important result of this study is the identification 
of an optical signature for the presence of a starburst in these 
galaxies, in spite of their large dust attenuation. While con- 
tinuous star formation follows a line along decreasing Balmer 
EW and increasing D4000, a burst superimposed on this pop- 
ulation produces a loop that first decrease D4000 and then in- 
crease the Balmer line EW. However, after about 1 .5 Gyr, the 
memory of the burst is lost and the galaxy behaves like others 
which did not experience a starburst. As result, we are limited 
to study only the averaged past star formation history for look- 
back times shorter than 1 .5 Gyr. The bursts characteristics were 
derived from probability distribution functions (PDF) using a 
bayesian statistics as in Kauffmann et al. (2003). 

The median ratio of present over averaged SFR, the so- 
called Scalo parameter, for the distant LIRGs is SFR/<SFR>= 
4±1 (we used the H8 line for which a larger sample of galax- 
ies is available and the PDFs present a sharper peak), which 
indicates that these galaxies are experiencing an atypically in- 
tense phase of star formation in their lifetime. A median SFR 
of 52^3 M yr~' for the ongoing starbursts was derived from 
their MIR luminosities. Hence, their mean SFR averaged over 
their lifetime is <SFR>= 13+3 M yr _1 . Knowing the median 
stellar mass for LIRGs of equivalent luminosity and redshift 
range (from Franceschini et al. 2003), we derived an age for 
those distant LIRGs of ff or m = 5.5+1.5 Gyr, suggesting that 
they formed at Zf orm = 4.5 + 1.5. 

As for the bursts themselves, we computed a median du- 
ration of t b = x 10 8 years, during which the galaxies 
produced 10±9% (the error bar includes 68 % of the galaxy 



sample) of their stellar mass. This corresponds to a mass of 
molecular gas of about 7xl0 9 M Q which is consistent that ob- 
served in local LIRGs and ULIRGs (see Sanders & Mirabel 
1996). 

We note that all simulations produced in this paper assume 
the same fixed IMF for both the underlying star formation 
and the burst of star formation. Some evidence that the for- 
mation of low-mass may be less efficient in the environment 
of active star formation in the solar neighborhood were sug- 
gested in the past (Larson 1986, Scalo 1986, Maeder 1993). 
A top-heavy IMF could also account for the enhanced ratio of 
light elements to iron in massive early-type galaxies (Worthey, 
Faber & Gonzalez 1992) and for the relative enrichment of oxy- 
gen to iron in the intra-cluster medium (Arnaud et al, 1992). 
In our study, the occurrence of a top-heavy IMF in the burst 
episode would only weakly influence the derived burst duration 
timescales, which are set by the spectral signature of massive 
A to F stars. However, a top-heavy IMF during the burst phase 
would imply a much lower contribution to the total galaxy mass 
by longer-lived, low-mass stars. 

Finally, we discussed the possibility that the distant LIRGs 
experienced a previous starburst prior to the ongoing one dur- 
ing the last 1 .5 Gyr. While most galaxies are not consistent with 
a merger scenario where two galaxies merge in several phases 
producing a series of bursts separated by a few ten millions 
years (Mihos & Hernquist 1996), the majority present more 
than 50 % chance to have experienced a previous burst in the 
last 1.5 Gyr, i.e. since z ~ 1. If these properties are typical of 
LIRGs between z = 1 and z = 0, then this suggests that the 
population of galaxies experiencing LIRG phases experienced 
on average 2 to 3 LIRG phases since z = 1 and up to 4 since 
their birth around z = 4-5, as also suggested by Hammer et al. 
(2005). This scenario is not consistent with the formation of 
distant LIRGs through the continuous star formation character- 
izing isolated spiral galaxies as has been argued independently 
based on their morphology. Instead, minor mergers, tidal inter- 
actions and gas accretion all remain equally plausible trigger- 
ing mechanisms for more than half of the distant LIRGs which 
do not harbor the morphology of major mergers. 
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